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SECTION  1 


INTRODUCTION 


For  many  purposes,  it  is  quite  useful  to  have  a simple  analytic 
model  for  estimating  SGEMP  effects.  Any  simple  model  will  undoubtedly  be 
somewhat  inaccurate  yet  such  a model  may  still  be  quite  valuable  in  pointing 
out  general  trends  and  in  promoting  a basic  understanding  of  SGEMP  phenomena. 
This  report  discusses  one  such  simple  model. 

This  simple  model  is  based  upon  a sphere  emitting  electrons  in  an 
azimuthal ly  symmetric  pattern  over  some  range  of  polar  angles.  By  assuming 
that  surface  charge  immediately  redistributes  over  the  sphere,  one  can  use 
the  continuity  equation  to  write  an  analytic  expression  for  the  surface 
replacement  current  in  terms  of  the  electron  emission  current  density.  These 
expressions  will  be  derived  and  the  results  then  compared  to  more  accurate 
numerical  calculations. 

These  analytic  expressions  for  replacement  current  will  then  be  used 
to  point  out  possible  difficulties  with  SGEMP  replacement  current  simulation 
techniques  where  only  a limited  number  of  drive  points  are  used  to  excite 
the  test  object.  It  is  shown  that  it  may  be  quite  difficult  to  obtain  the 
proper  skin  currents  on  the  emitting  surface  of  the  test  object  unless  the 
entire  surface  is  excited, as  in  the  case  of  SGEMP.  For  skin  currents  in 
the  shadowed  region,  however,  the  effect  of  discrete  drivers  is  not  nearly 
so  important. 


SECTION  2 


A SIMPLE  MODEL  FOR  ESTIMATING  REPLACEMENT  CURRENTS 


The  basic  geometry  considered  is  spherical,  as  indicated  in  Figure 
1.  Azimuthal  symmetry  with  respect  to  the  z-axis  is  assumed.  Some  radial 
emission  current  density,  J,  is  also  assumed;  i.e., 

J = J0(r,t)[U(ei)  - U(02)]  i , (1) 

/v 

where  U(0)  is  the  unit  step  function  and  r is  a unit  vector  in  the  radial 
direction. 


The  basis  for  this  simple  model  is  the  continuity  equation, 

||  + V • J » 0 , (2) 

which,  when  calculated  at  the  surface  of  the  sphere  becomes 

dd~  = Jotr=a,0)  + a sin0  w(sin6  K(0^)  1 

2 

where  a is  the  surface  charge  density  (coul/m  ) and  K is  the  surface  current 
density  (amps/m).  Note  that  a and  K are  just  the  basic  SGEMP  response 
parameters . 

In  general,  one  must  solve  both  Maxwell's  equations  and  the  equa- 
tions of  motion  of  the  emitted  electrons  to  solve  for  the  a and  K parameters 
in  Equation  3.  Several  simplifications  to  this  process  are  possible,  however. 
First  of  all,  one  can  specify  the  spatial  current  density,  J.  This  is 


► 
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2 
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Figure  1.  Basic  geometry. 

equivalent  to  the  assumption  that  electron  trajectories  are  fixed  and  will 
not  be  influenced  by  the  resultant  electromagnetic  fields.  Further 
simplification  is  possible  if  one  assumes  the  problem  is  quasi-static.  One 
can  then  solve  Poisson's  equation  for  the  electric  field,  given  the  boundary 
conditions  and  electron  positions  in  space.  The  surface  charge  density,  a, 
can  then  be  calculated  from  the  electric  field  normal  to  conducting  surfaces 
and  Equation  3 then  gives  the  surface  current  density,  K. 

Note  that  an  accurate  quasi-static  calculation  includes  the  ef- 
fects of  electrons  that  have  been  emitted  from  the  test  object  but  have  not 
yet  had  time  to  move  very  far  from  the  object.  The  presence  of  these  nearby 
electrons  tends  to  hold  down  surface  currents  since  they  repel  any  electrons 
flowing  on  the  surface  to  replace  the  emitted  charge.  [One  can  also  describe 
this  process  by  saying  that  the  image  charge  has  not  yet  moved  away  from  the 
emission  surface. ] 

A worst  case  estimate  of  skin  currents  can  thus  be  obtained  by 
assuming  that  all  emitted  electrons  immediately  move  to  infinity.  In  this 
case,  the  resultant  skin  currents  are  indeed  the  so-called  "replacement 
currents"  since  charge  will  immediately  flow  to  replace  emitted  electrons. 
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Under  these  conditions,  the  surface  charge  density  is  just  that 
which  would  exist  if  the  object  were  DC  charged.  For  the  conducting  sphere, 
this  charge  density  is  uniform;  i.e., 

0(0)  = , (4) 

4ira 

where  Q is  the  total  charge  on  the  sphere.  This  equation  can  be  rewritten 
as 


do  _ *E  J0Ae 
dt  ' " \ ' \ ’ 


(5) 


where  I_  is 
2 E 

(4ira  ) and  A 


the  total  emission  current,  Aj,  is  the  total  area  of  the  sphere 
is  area  of  that  part  of  the  sphere  that  is  emitting  electrons. 


If  Equations  3 and  5 are  combined,  remembering  that  negative 
charge  (i.e.,  electrons)  is  emitted,  it  can  be  shown  that 

|j0|Ae  2 

1,(9)  —-5-  2na^  (1  - cos6)  for  0 s 6 < 0j  , (6) 


and 


where 


A A 

V6>  ‘ • ^Ae  lJ0l  * M’  - a;)2™2^6!  - cos6) 


for  e s e < e2  , 


(7) 


i‘J0*Ae  2 

rK(0)  = 2ira  (1  + cos6)  for  02  < 6 < tt 


(8) 


Ax  = 2Tta“(l  - cosex)  , 


Ag  = 2na  (cosO^  - cos62)  , 


(9) 

(10) 


and 
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IK(0J  = 2ita  sin6  K(0)  . 


(ID 


Note  that  1^(0)  is  just  the  total  replacement  current  as  a function  of 
angle  0. 


Then 


This  general  case  can  be  simplified  if  it  is  assumed  that  6^  = 0. 


IK(9)  = !J0l(1  ' A^)2lTa2(1  "cos9)  for  0i  6 <02  , (12) 


and 


lJ0 lAe  2 

!k(8)  = — ^ 2Tra^(l  + cos0)  for  02  S 0 < Tr  . 


(13) 


These  two  eqt ations  indicate  several  very  interesting  facts. 

First  of  all,  the  peak  SKin  current  will  occur  at  0 = 0^,  the  edge  separat- 
ing the  emitting  and  non-emitting  regions  of  the  surface.  The  peak  current 
is  then 


h = lJol(i  - r)  27ra2(i  - cos9 2)  * 

PEAK  ’ \ ' 

However, 

A = 2Tra2(l  - cos9_)  , 
e 2 


(14) 


(15) 


so  that 


(16) 
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where  1^  is  the  total  emitted  current. 

Equation  16  says  that  the  peak  replacement  current,  for  a given 
emission  current,  just  depends  upon  the  fraction  of  the  surface  area  that  is 
emitting  electrons!  Furthermore,  when  the  whole  surface  is  emitting,  replace- 
ment currents  go  to  zero.  [This  result  is  quite  reasonable  since  the 
problem  becomes  spherically  symmetric  when  the  entire  surface  is  emitting 
and  it  is  well  known  that  for  spherical  symmetry  only  radial  electric  fields 
are  created;  i.e.,  no  magnetic  fields  or  skin  currents  are  produced.]  It 
is  also  obvious  that  the  maximum  value  of  replacement  current, for  a given 
emission  current,  is  achieved  when  the  emission  area  is  small  (i.e.,  when  the 
emission  current  is  emitted  along  a filament  along  the  z-axis). 


For  most  SGEMP  cases,  one  would  expect  about  half  the  external 
surface  to  be  emitting  electrons  (assuming  the  body  is  not  highly  convoluted) . 
Therefore,  the  peak  replacement  current  can  be  very  roughly  estimated  by 


[ < -y- 

kpeak~  2 


One  should  again  be  reminded  that  this  is  a worst  case  estimate 
and  that  actual  replacement  currents  will  depend  on  electron  emission  spectra 
and  the  average  speed  that  the  emitted  electron  cloud  moves  away  from  the 
emission  surface. 


SECTION  3 


COMPARISON  WITH  MORE  ACCURATE  CALCULATIONS 


Skin  currents  on  a 1 meter  radius  sphere,  as  calculated  by  the 
simple  scheme  described  here  and  as  obtained  from  the  LFLUX  code1,  are 
compared  in  Figure  2.  The  prescribed  emission  current  had  a sin2(“ir) 
time  history  with  T = 10  ns  and  a peak  emission  current  density  of  1.0 
amps/m  . The  LFLUX  calculation  assumed  all  electrons  were  emitted  radially 
from  half  the  spherical  surface  and  that  the  electron  soectrum  was  mono- 

g 

energetic,  corresponding  to  a velocity  of  1.0  x 10  m/sec.  Note  that  the 
LFLUX  code  solves  Maxwell's  equations  numerically,  given  a prescribed 
current  density. 

A comparison  of  the  currents  predicted  by  Equation  12  and  by  LFLUX 
(assumed  here  to  be  the  baseline  for  evaluating  the  accuracy  of  our  simple 
model)  shows  that  the  simple  model  does  give  : eplacement  currents  several 
times  larger  than  more  accurate  calculations.  The  reason  for  this  error 
is  the  fact  that  emitted  electrons  are  assumed  to  immediately  move  to 
infinity,  thus  allowing  complete  charge  redistribution  immediately  after 
electron  emission.  On  the  other  hand,  the  LFLUX  calculation  includes  the 
fact  that  the  electrons  take  some  finite  time  to  move  away  from  the  sphere, 
thus  giving  a lower  peak  current  and  a current  pulse  stretched  out  in  time. 


1 Stettner,  R.,  and  0.  Higgins,  X-ray  Induced  Currents  on  the  Surface  of 
a Metallic  Sphere,  Mission  Research  Corporation,  MRC-N-111,  DNA  3612T, 

April  1975. 
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Time  (ns) 


Figure  2.  A comparison  of  replacement  current 
densities  as  calculated  from  Equation 
12  and  from  the  LFLUX  code. 
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Note  that  the  total  replacement  charge  (i.e. , the  time  integral 
of  the  replacement  current)  will  be  the  same  for  LFLUX  and  the  simple 
analytic  model  used  here  as  long  as  all  of  the  emitted  electrons  eventually 
escape  the  sphere  and  move  to  infinity.  Thus,  the  simple  model  developed 
here  will  always  give  a replacement  current  with  too  high  of  a peak  value 
and  too  short  of  a pulse  width.  The  product  of  pulse  height  and  pulse 
width  will,  however,  be  roughly  the  same  as  that  of  a more  accurate  calcula- 
tion (at  least  for  triangle-type  time  histories  where  this  product  is  just 
proportional  to  the  total  replacement  charge) . In  mathematical  terms 


Tr  = I T “0 
t t 0*0  y 


(18) 


where 

It  is  the  peak  current  as  given  by  Equation  12  or  13 

Tt  is  the  pulse  width  of  the  emission  current 
(e.g.,  FWHM) 

IQ  is  the  actual  peak  replacement  current  flowing 
past  some  location 

Tq  is  the  actual  pulse  width  of  the  replacement 
current  pulse 

and 

Q is  the  total  replacement  charge  that  flows  past 
the  specified  location. 


Note  that  for  0 = 90°,  Figure  2 gives 
I = 3.14  amps 
Tt  = 5 ns 
Iq  = 1.31  amps 
Tq  = 9.8  ns 
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so  that  Equation  18  is  valid  to  within  about  20  percent. 


Another  way  of  looking  at  the  simple  analytic  model  being  con- 
sidered here  is  to  compare  it  to  a circuit  model  that  treats  the  space  out- 
side the  sphere  as  a single  capacitor  with  a parallel  current  generator 
representing  the  electron  emission  current.  It  is  known2  that  such  a model 
gives  replacement  currents  from  two  to  ten  times  larger  than  a more  accurate 
circuit  which  uses  numerous  radial  capacitors  and  takes  the  finite  electron 
velocity  into  account. 


To  correct  for  this  effect,  one  can  also  calculate  an  "effective" 
capacitance.  If  at  some  time  t,  the  "center"  of  the  emitted  electron  cloud 
is  some  distanc  b from  the  center  of  a sphere  of  radius  a,  the  effective 
capacitance  is  just 


4tte  ab 

p = z 

eff  b - a • 


(19) 


[Cfiff  is  just  the  capacitance  of  concentric  spheres 


If  the  electrons  had  moved  to  infinity,  however,  the  capacitance 
is  given  by 

c«  = 47re0a  ' (20) 

The  amount  of  charge  redistribution  is  then  just  proportional  to 

the  ratio 


(21) 


Note  that  this  ratio  equals  .5  when  b = 2a;  i.e.,  half  of  the  total  charge 


2 Wenaas,  E.  P. , "Lumped -Element  Modeling  of  Satellite  SGEMP  Excitation," 
IEEE  Trans.  Nucl.  Sci.,  Vol.  NS-21,  December  1974. 
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replacement  will  have  occurred  when  the  "average"  electron  has  moved  one 
spherical  radius  away  from  the  emission  surface. 


One  can  thus  roughly  estimate  the  actual  replacement  current  pulse 
width  T0  (FWHM)  by 


(22) 


where  v^  is  the  average  electron  velocity.  [This  formula  is  valid  only 
for  emission  time  histories  with  pulse  widths  less  than  T^.] 


Equation  18  can  be  combined  with  Equation  22  to  give  an  estimate 
of  the  actual  peak  replacement  current. 


I 


0 


(23) 


Equations  22  and  23  agree  fairly  well  with  the  data  shown  in  Figure  2. 
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SECTION  4 


SIMULATION  CONSIDERATIONS 


Now  let  us  use  Equations  6,  7,  and  8 to  look  at  the  distribution 
of  surface  replacement  currents  around  the  sphere.  Figure  3 shows  the 
normalized  replacement  current,  l^/2n,  as  a function  of  polar  angle  0 on  a 
1 meter  radius  sphere  with  an  emission  current  density  JQ  of  1 amp/m2. 

Four  different  cases  are  shown:  Case  1 is  for  emission  in  the  range 
0°  s 0 < 90°,  Case  2 for  0°  < 0 < 30°,  Case  3 for  30°  < 0 <60°,  and 
Case  4 for  60°  < 0 < 90° . 

Note  that  the  current  for  Case  1 is  just  the  summation  of  the  cur 
rents  for  Cases  2,  3,  and  4,  as  expected  since  the  superposition  of  these 
emission  regions  is  just  the  emission  region  of  Case  1.  Note,  however, 
tnat  there  is  significant  cancellation  of  currents  when  Cases  2,  3,  and  4 
are  summed  in  the  range  0°  < 8 < 75°.  This  is  due  to  the  fact  that,  for 
example,  emission  from  60°  to  90°  only  gives  negative  skin  currents  for 
0“  ~ 9 - 60°  while  half  sphere  emission  gives  positive  currents  in  this 
region.  In  Case  4,  electrons  will  flow  from  the  front  of  the  sphere  to 
replace  those  emitted,  while  in  Case  1,  electron  emission  from  th3  front  of 
the  sphere  will  dominate  this  effect  thus  giving  rise  to  currents  of  the 
opposite  sign. 

This  cancellation  effect  creates  some  interesting  questions 
regarding  various  simulation  techniques  that  may  be  used  to  test  satellites 
for  SGEMP  response.  For  example,  current  sources  may  be  attached  to  a 
satellite  at  a number  of  locations  with  thin  wires  to  simulate  electron 
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emission.  Since  only  a finite  number  of  current  generators  would  be  used 
to  simulate  electron  emission  over  an  entire  exposed  surface,  some  ques- 
tion occurs  as  to  the  accuracy  of  the  simulation  technique;  i.e.,  how  many 
drive  points  are  required  and  what  timing  requirements  are  necessary  to 
reproduce  the  threat  situation? 


For  example,  instead  of  individual  wires,  let  us  assume  we  have 
azimuthully  symmetric  current  sources  and  that  we  wish  to  simulate  half- 
sphere emission  (0°  < 0 £ 90°)  by  two  drivers,  one  extending  over 
0°  s 0 < 30°  and  the  other  in  the  range  60°  < 0 s 90°.  For  a 1 meter  radius 
sphere,  the  surface  area  for  9^  £ 0 5 0^  is 

A(01  < 0 < 02)  = 2ir(cos01  - cos92)  . (24) 

Therefore 

A(0°  < 0 < 30°)  = .84  m2  = AJ 

A(60°<  0 < 90°)  * 3.14  m2  * A2  (25) 

A(0°  < 0 < 90°)  = 6.28  m2  = A3 

Let  us  ssume  that  the  emitted  currents  are  scaled  by  the  ratio 


A1  + A2 


= 1.58  , 


(26) 


so  that  proper  total  emitted  current  is  obtained.  One  can  then  multiply 
the  appropriate  curves  in  Figure  3 by  this  ratio  and  add  the  results  to 
obtain  the  simulated  skin  current  distribution.  The  result  of  this 
operation  is  shown  in  Figure  4. 


An  examination  of  Figure  4 reveals  several  interesting  facts. 
First  of  all,  the  "simulated"  skin  current  distribution  does  a good  job 
of  reproducing  the  actual  skin  current  distribution  on  the  back  half  of 
the  sphere;  i.e.,  for  90° * 9 S 180°.  On  the  front  surface  (0°  ^ 0 < 90°), 
however,  simulation  is  quite  poor,  with  the  "simulated"  current  exceeding 
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Polar  Angle  0 (degrees) 

Figure  4.  Comparison  of  replacement  current  distributions  for  Case  1 (half  sphere 
emission)  and  "simulated1’  excitation. 


the  actual  value  for  0°  - 0 £ 40°  and  dropping  below  the  actual  value  for 
40°  - 0 £ 90°.  At  6 - 60°  the  "simulated"  current  even  has  the  wrong  sign! 

The  importance  of  such  deviations  will  depend  upon  the  system  being 
tested.  If  an  aperture  or  antenna  was  located  at  the  0 = 60°  location,  this 
simulation  technique  might  give  considerable  error  in  evaluating  the 
amount  of  coupling  to  the  inside  of  a system.  If  the  important  points-of- 
entry  were  on  the  back  of  the  sphere,  however,  relatively  good  test  data 
should  be  attainable,  however. 

The  effect  of  having  a local  minima  or  null  on  the  surface  between 
two  drive  points  is  expected  to  be  a real  physical  phenomena  rather  than 
one  dependent  upon  the  simple  model  used  here  to  demonstrate  the  effect. 

Since  any  two  drive  points  will  be  simulating  electron  emission  from  the 
surface,  it  is  quite  reasonable  to  expect  some  point  on  the  surface  between 
the  two  sources  that  is  equally  influenced  by  each,  resulting  in  a null! 

One  must  thus  be  quite  careful  in  applying  SGEMP  simulation  techniques  to 
areas  that  are  actually  emitting  electrons. 
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SECTION  5 


4 


CONCLUSIONS  AND  SUMMARY 


In  this  report,  a simple  analytic  model  for  estimating  SGEMP 
replacement  currents,  given  the  emission  current,  is  derived  and  compared 
with  a more  accurate  numerical  calculation.  The  simple  model  is  shown  to 
overestimate  the  peak  replacement  current  magnitude  and  underestimate  the 
pulse  width.  Despite  these  limitations,  this  simple  model  is  still  quite 
useful  for  making  rough  estimates  and  for  promoting  a general  understanding 
of  replacement  current  phenomena. 

This  simple  model  was  then  used  to  demonstrate  a potential  problem 
with  SGEMP  simulation  techniques  that  drive  only  discrete  locations  rather 
than  an  entire  surface.  Such  drive  techniques  are  shown  to  give  good 
simulation  of  replacement  currents  on  the  rear  (undriven)  part  of  the  test 
object,  but  unrealistic  local  skin  current  minima  are  shown  to  be  produced 
at  some  point  between  injection  points. 
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ATTN:  Conrad  L.  Longmire 
ATTN:  Daniel  F.  Higgins 
ATTN:  Roger  Stetlner 
ATTN:  William  C.  Hart 
5 cy  ATTN:  Doc.  Con. 

Mission  Research  Corporation 
ATTN:  David  E.  Merewether 
ATTN:  Larry  D.  Scott 

Mission  Research  Corporation-San  Diego 
ATTN:  V.  A.  J.  Van  Lint 

The  Mitre  Corporation 

ATTN:  Theodore  Jarvis 

Motorola,  Inc. 

Government  Electronics  Division 

ATTN:  James  R.  Black,  M.S.  A-112 
ATTN;  A.  J.  Kordalcwski,  Tech.  Info.  Ctr. 

Northrop  Corporation 
Electronic  Division 

ATTN:  John  M.  Reynolds 
ATTN:  Vincent  R.  DeMartino 
ATTN:  Boyce  T.  Ahlport 
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Northrop  Corporation 
ATTN:  Library 
ATTN:  Orlie  L.  Curtis,  Jr. 
ATTN:  David  N.  Pocock 


Sanders  Associates,  Inc. 

ATTN:  R.  G.  Dcspathy,  Sr.,  PE,  1-6270 
ATTN:  Moe  L.  Altel,  NCA,  1-3236 
ATTN:  James  L.  Burrows 


Northrop  Corporation 

Electronic  Division 

ATTN:  Joseph  D.  Russo 

Physics  International  Company 

ATTN:  Doc.  Con.  for  Philip  W.  Spence 
ATTN:  Doc.  Con.  for  John  tl.  Huntington 
ATTN:  Doc.  Con.  for  lan  D.  Smith 
ATTN:  Doc.  Con.  for  Charles  H.  Stallings 

Pulsar  Associates,  Inc. 

ATTN:  Carleton  H.  Jones,  Jr. 

RAD  Associates 

ATTN:  William  J . Karzns 
ATTN:  S.  Clay  Rogers 
ATTN:  Richard  R.  Schaefer 
ATTN:  Leonard  Schlesslngcr 
ATTN:  William  R.  Graham,  Jr 

The  Rand  Corporation 
ATTN:  Cullen  Crain 

Raytheon  Company 

ATTN:  Gajanan  fl.  Joshi,  Radar  Sys.  Lab. 

Raytheon  Company 

ATTN:  Harold  L.  Flo  sc her 
ATTN:  James  R.  Weckback 

RCA  Corporation 

Government  A Commercial  Systems 
ATTN:  George  J.  Brucker 

RCA  Corporation 

ATTN:  K.  H.  Zalntngcr 

RCA  Corporation 

Government  & Commercial  Systems 
ATTN:  Andrew  L.  Warren 
ATTN:  Eleanor  K.  Daly 

RCA  Corporation 

ATTN:  E . Van  Keuren,  13-5-2 

Research  Triangle  Institute 

ATTN:  Sec.  Officer  tor  Eng.  Dlv  , Mayrant  Simons,  Jr. 

Rockwell  International  Corporation 
ATTN;  Donald  J.  Stevens,  FA-70 
ATTN:  L.  H.  Pinson,  FB-41 
ATTN:  James  E.  Bell,  HA-10 
ATTN:  J.Spetz 
ATTN:  N.J.Rudie,  FA-53 
ATTN:  George  C.  Messenger,  FB-01 


Science  Applications,  Inc. 

ATTN:  William  L.  Chadse.v 

Science  AppI ications,  Inc. 

ATTN:  Noel  R.  Byrn 

Science  Applications,  Inc. 

ATTN:  J.  Roger  Hill 

Science  Applications,  Inc. 

ATTN;  Charles  Stevens 

Sidney  Frankel  £ Associates 
ATTN:  Sidney  Frankel 

Simulation  Physics,  Inc. 

ATTN:  Roger  G.  Little 

The  Singer  Company 

ATTN;  Irwin  Goldman,  Eng.  Management 

Sperry  Flight  Systems  Division 
ATTN:  D.  Andrew  Schow 

Sperry  Rand  Corporation 
Univac  Division 

ATTN:  James  A.  mda,  M.S.  41T2S 

Sperry  Rand  Corporation 
Sperry  Division 

ATTN:  Paul  Marraffino 
ATTN:  Charles  L.  Craig.  EV 

Stanford  Research  Institute 

ATTN:  Robert  A.  Armistead 
ATTN:  Me!  Bernstein 
ATTN:  Setsuo  Dalrlkl 
ATTN:  Philip  J.  Dolan 
ATTN:  Arthur  Lee  Whitson 

Stanford  Research  Institute 

ATTN:  MacPherson  Morgan 

Sundslrand  Corporation 

ATTN:  Curtis  B.  White 

Systems,  Science  i Software 
ATTN:  David  A,  Dleskan 

Systems,  Science  i Software,  Inc. 

ATTN:  Ira  Katz 
ATTN:  Andrew  R.  Wilson 

Systron-Donner  Corporation 
ATTN:  Harold  D.  Morris 


Rockwell  International  Corporation 
ATTN;  John  F.  Roberts 

Rockwell  International  Corporation 
ATTN:  T.  B.  Yates 


Texas  Instruments,  Inc. 

ATTN:  Donald  J.  Manus,  M.S.  72 

Texas  Tech  University 

ATTN:  Travis  L.  Simpson 
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TRW  Systems  Group  United  Technologies  Corporation 

ATTN:  William  H.  Robinette,  Jr.  Norden  Division 

I ATTN:  Tech.  Info.  Ctr„  S-1930  ATTN:  Conrad  Corda 

ATTN:  Donald  Jortncr 

ATTN:  Richard  H.  Kingsiand,  Rl-2154  Westinghouse  Electric  Corporation 

ATTN:  Robert  M.  Webb,  111-1150  ATTN:  Henry  P.  Kulapaoa,  M.S.  3525 

ATTN:  Jerry  1.  Lubeli 

ATTN:  Aaron  H.  Narevsky,  Rl-2144 

ATTN:  Paul  Molmud,  R1-U96 

ATTN:  Philip  R.  Carder,  Rl-1028 

ATTN:  Lillian  D.  Singletary,  Rl-1070 

TRW  Systems  Group 

San  Bernardino  Operations 
ATTN:  J.  M.  Gorman 
ATTN:  John  E.  OahnKe 
ATTN:  Earl  W.  Allen 

TRW  Systems  Group 

ATTN:  Donald  W.  Pugstey 
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